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ABSTRACT 

Gamma-ray bursts (GRBs) that emit photons at GeV energies form a small but significant population 
of GRBs. However, the number of GRBs whose GeV-emitting period is simultaneously observed in 
X-rays remains small. We report 7-ray observations of GRB 110625 A using Fermi's Large Area 
Telescope (LAT) in the energy range 100 MeV-20 GeV. Gamma-ray emission at these energies was 
clearly detected using data taken between 180 s and 580 s after the burst, an epoch after the prompt 
emission phase. The GeV light curve differs from a simple power-law decay, and probably consists 
of two emission periods. Simultaneous Swift/XRT observations did not show flaring behaviors as 
in the case of GRB 100728A. We discuss the possibility that the GeV emission is the synchrotron 
self-Compton radiation of underlying ultraviolet flares. 
Keywords: gamma rays: bursts — gamma rays: observations 



1. INTRODUCTION 

Since the launch of the Fermi satellite in 2008, more 
than 20 7-ray bursts (GRBs) have been detected above 
~100 MeV by the Large Area Telescope (LAT) aboard 
the satellite (|Abdo et al.ll2009allbHc]fd1 '). Long-lived MeV- 
GeV emission of GRBs, first detected in the EGRET era, 
is now a common feature of LAT-detected GRBs. The 
nature of such temporally extended emission beyond the 
prompt GRB phase is not well understood. A widely- 
discussed radiation mechanism is sync hrotron emission 
from external shocked e l ectron s (e.g.. [Zou et al. 120091 : 
iKumar fc Barniol Duranl 120091 : iGhisellini et all I2010D . 
but inverse-Compton scattering off flare photons or late- 
time activ ities of the central engine are among a lternative 
scenarios (jAbdo et al.ll2011tlZhanglt^[2lFll . 

Simultaneous observations at other wavelengths of 
such extended MeV-GeV emission from GRBs are 
crucial to disentangle various emission models. By 
May 2011, only two LAT-detected GRBs have been 
simultaneously observed by Swift's X-ray telescope 
(XRT) during its GeV-emitting period: GRB 090510 
and GRB 100728A. GRB 090510 remains the only 
short GRB detected by LAT, and its GeV emission 
can be interpreted as, e.g., synchrotron radiat i on of 
the forward shock el ectrons (iDe Pasauale et ell 120101 : 
Ghirlanda et~ap I2010L but see iGao et all (|2009h and 
Piran fc Nakarl (|2010l ) for an opposing viewpoint). In 



the case of GRB 100728A, X-ray flares and correspond- 
ing GeV emission were detected up to ~1 ks after 
the b urst, suggesting their common origin ([Abdo et al.l 
1201 If ) . However, the number of GRBs whose GeV emis- 
sion is simultaneously covered in X-rays remains low. 

In this paper, we report another such case: 
GRB 110625 A, that was detected by Fermi/LAT and 
Swift/XRT simultaneously for several hundred seconds. 
Errors are reported at ler, unless otherwise specified. 
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2. GRB 110625A 

At 21:08:28 UT on 2011 June 25, the Burst 
Alert Telescope (BAT) aboard Swift triggered on 
GRB 110625A (jPage et al.l l2011al ). The refined BAT 
position was R.A. = 19 h 07 m 00!3, Dec. = +06°45T7'.'8 
(J2000) with an uncertaint y of 1.'3 (90% containment 
radius; iPalmer et all l2011h . The BAT light curve 
showed a multiple-peaked structure lasting from ~ — 12 s 
to ~ 18 s with a tail extending up to ^150 s with 
respect to BAT trigger time. T 90 (referred to the 
time interval between the instants at which 5% and 
95% of the total integral emissio n is detected in th e 
15-350 keV ba nd) was 44.5 ± 10. 1 s (IPalmer et al.l[20Tl . 
Konus-Wind (|Golenetskii et~aT1 l2011lh Fermi/GBM. 
INTEGRAL S PI-ACS (IGruber et all IMtl . and 
Suzaku/WAM (|Mizuno et al.l I2011D also triggered 
on GRB 110625A. The Konus-Wind team reported 
that the 20 keV to 10 MeV time-averaged spec- 
trum from to 58.88 s after the Konus-Wind trig- 
ger time is best fitted by the Band function with 
a = -1.05 ± 0.08, (3 = -2.7±g;§, and E p = 190±^ keV; 
emission was seen up to ~8 keV. The burst fluenc e is 
(6.1 ± 0.6) x 10" 5 erg cm" 2 (jGolenetskii et al.|[20TTh . 

Fermi/GBM triggered on GRB 110625A at 21:08:18.24 
UT (To) on 2011 June 25. The angle of the GRB po- 
sition is 88° from the LAT boresight at To- The au- 
tonomous rapid repoint maneuver repointed the LAT 
such that GRB 110625A was put in the field-of-view 
(FoV) of LAT from - T + 100 to T + 600. However, 
due to the poorly measured position derived on board 
by GBM (off by ~68°), The burst position was placed 
at the outs kirt of the FoV of Fe rmi/LAT, diminishing its 
sensitivity (jGruber et al.ll2011l ). 

Swift/XRT began data-taking of the burst at « To + 
150 s and found the X-ray afterglow source at the po- 
sition R.A = 19 h 06 m 55!85, Dec. = +6°42T9'.'2 (J2000; 
IPagdl20i"ll ) with an error circle of radius 2'.'1 (90% confi- 
dence level). This position is used in analyses presented 
in this paper. The flux faded initially as a power law with 
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an index of ot\ — 1.14 ± 0.04 and then as a-i = 2.3 jlo'4, 
after a break at 17^ ks after th e BAT trigger (90% 
confidence level; iPage et alJfeOllbl ). 

The XRT observations of GRB 110625 A started in 
Window Timing (WT) mode that lasted for «90 s. Then 
data were taken in the Photon counting (PC) mode 
since ~ T + 250 s. The WT(PC) spectrum can be 
fit by an absorbed power-law model with photon inde x 
2.5 ± 0.4(1.8 ± 0.4) (90% confidence level; [Pa^[20Tl . 
suggesting a marginal hardening between early and late 
times. We tried the absorbed blackbody model as well, 
finding no improvement on the fit on both WT- and PC- 
mode data. Using PC-mode data taken from To + 255 s 
to To + 580 s (to match LAT observations), we fit the 
time-averaged 0.3-10 keV energy spectrum with an ab- 
sorbed power- law model. A photon index of 1.44+o'22 
and a hydrogen column density of 4. 53^59 X 10 22 cm -2 
were obtained. In our analysis, we used C-statistic dur- 
ing the fitting process. Piled-up effects are present in the 
PC-mode data concerned and were removed by ignoring 
the circular region with radius 7.'2. 

Interestingly there is no clear early steep decay phase, 
a common feature in Swift XRT afterglows that has been 
widely interpreted as the high latitude emission from the 
prompt (Fenimore et al. 1996; Kumar & Panaitescu 
2000). The absence of a sharp decline phase may be 
just due to the fact that XRT observations started only 
at « T + 150 s, i.e., the afterglow started earlier. Such a 
possibility provides us an additional argument to classify 
the LAT emission as post prompt emission. 

In the optical and infrared bands, the only object sus- 
pe cted of being var i able i n the XRT FoV was reported 
in iGorosabel et al.l (|2011l ). However, from an archival 
Canada- France-Hawaii Telescope r-band image (taken on 
2010 May 8), the object was detected with r = 23.9±0.1 
and not detected with the g-band filter with a 5-sigma de- 
tection l imit of 26.2, consis tent with the GROND obser- 
vations (Filga s et al.l 1201 ll ). We therefore conclude that 
it is unrelated to GRB 110625A. 

3. THE FERMI/LAT DATA ANALYSIS AND RESULTS 

We analyzed the LAT data that are available at the 
Fermi Science Support CenteiQ. The Fermi Science Tools 
v9r23pl package was used to reduce and analyze the 
data. We selected photons of energies between 100 MeV 
and 20 GeV. To reduce the contamination from Earth 
albedo 7-rays, we excluded events with zenith angles 
greater than 100°. 

We selected photons from a region-of-interest (ROI) of 
a 10°-radius circular region centered on the XRT position 
of GRB 110625A and plot them in Fig. [TJ Here we use 
"P7TRANSIENT" data so as to increase photon statis- 
tics. We selected the time span during which the inclina- 
tion angle of the GRB position is less than 66° to make 
sure that the GRB position is well within the LAT FoV, 
corresponding to T +180 s to T o +580 s. The count rate is 
normalized by the varying exposure, that was computed 
by using the best-fit spectral index T 7 = 2.7 obtained in 
gtlike (described below). The background level of a time 

bin i was calculated by Pbg,i = ~%2f=i bj, where bj repre- 
sents the probability that a particular photon j in time 

1 http: //fermi .gsfc .nasa.gov/ssc/ 



bin i comes from the Galactic diffuse emission, i.e., from 
the background, and Ni the number o f photons in the 
ROI. A similar approach was used by I Guillemot et al.l 
1)20111 ). This probability is calculated using gtsrcprob 
based on the best likelihood fit betwee n Tp + 180 s and 
T + 580 s, according to the method of lKerrl (|20111 ). 
The weighted photon flux from GRB 110625 A was cal- 



culated by i"GRB,i = w j> wnere Wj = 1 — bj rep- 

resents the probability that the photon originates from 
the GRB, normalized by the exposure. Thus, the total 
number of counts, Ni, is scaled by the factor PGRB,i/Ai 
to give rise to the weighted photon counts. The er- 
ror bars were scaled by the same ratio. The back- 
ground level for the weighted photon flux is estimated 
by Bi — 53 i=i w j x for each bin and was averaged 
over the whole period. It appears that the weighted light 
curve differs from a simple power law decay. To quantify 
this, we fit the weighted light curve with a simple power 
law between T + 180 s and T + 580 s: f = fo + at,t mea- 
sured in seconds, and found the best fit parameters to be 
/o = 5.2xl0~ 5 cm- 2 s- 1 and a = -4.3xl0~ 8 cm- 2 s- 2 , 
giving x 2 /d.o.f.= 16.79/8. Thus the light curve deviates 
from a simple power law decay at the level of 96.77%. 

The ratios Pgrb,i/ Pbg,i are shown in Fig. Q] as well. 
One can see that during T + 260 s to T + 340 s and T + 
460 s to T + 500 s, the contribution from GRB 110625 A 
is higher than the background, since Pgrb,[/Pbg,i > 1- 
As an alternative algorithm, we have also assigned each 
photon to either GRB 110625 A or background depend- 
ing on whether Wj > 0.5 (GRB 110625A) or wj < 0.5 
(background) and calculated the ratio of the number of 
photons from GRB 110625A and those from background. 
This method gives consistent results. 

The existence of two emi ssion epochs from 
GRB 110625A was fi rst noted in [Tarn fc Kond (|2011h 
and characterized by iGruber et al. ()201lD . Fitting the 
light curve with two Gaussian profiles on top of a baseline 
emission, the peaks and widths of the first and second 
emission epochs are given by t pi = To + 285.8 ± 10.0 s 
and Ai = 64.9 ± 22.3 s, and by t P2 = T Q + 483.8 ± 15.7 s 
and A2 = 59.0 ± 31.0 s, respectively. Therefore, the two 
emission epochs are roughly T + 220.9 s to T + 350.7 s 
and T + 424.8 s to T + 542.8 s. 

We then performed unbinned maximum-likelihood 
analyses (gtlike) of a 15°-ROI centered at the XRT 
position to characterize the detection significance and 
spectrum of the >100 MeV 7-rays from GRB 110625 A. 
Events that are classified as "P7SOURCE" and arrived 
between T o + 180 s and T +580 s are used. We subtracted 
the background contribution by including the Galac- 
tic diffuse model (gal_2yearp7v6_v0.fits). Contributions 
from other sources (including the isotropic background) 
are negligible and not included. The spectral index of the 
GRB and normalization values of all three components 
were set free in the likelihood analysis. The likelihood fit 
returned a best-fit spectral index of T-y = 2.7 ± 0.3 and 
a test- statistic (TS) value (jMattox et al.1 119961 ) of 52.9, 
corresponding to a detection significance of ^7a. Using 
gtfindsrc, we found the best-fit LAT position to be right 
ascension (J2000) = 286?51 and declination (J2000) = 
6?86 with statistical uncertainty of 0?44, which is con- 
sistent with the XRT position. The systematic uncer- 
tainty should be >0?6 based on GRB 080825C which 



Fermi LAT observations of GRB 11 0625 A 



3 



Table 1 

Spectral properties of the GeV emission during different periods. 



i - To (s) 


Photon Index 


Tii 


1: 180-350 


-2.6±0.3 


42.8 


II: 350-425 







III: 425-580 


-2.7±0.6 


8.5 


180-580 


-2.7±0.3 


52.9 



also occurred at large inclination angle at the GRB on- 
set (|Abdo et al.H2009bl) . 

Performing unbinned maximum-likelihood analyses 
(gtlike) of the same 15°-ROI for different time bins, 
we produced the background-subtracted light curve, as 
shown in Figure [5] The first and third data point repre- 
sents the periods (I) T + 180 s to T + 350 s and (III) 
T + 420 s to T + 580 s, respectively. These two peri- 
ods, (I) and (II), were chosen to cover the two emission 
epochs characterized above, as well as the short periods 
T + 180 s to T + 220.9 s and T + 542.8 s to T + 580 s, 
respectively. Limited photon statistics do not allow any 
meaningful analysis on these ~40 s periods. The period 
in-between (II) during which no emission was detected is 
plotted as a 90% confidence-level upper limit. The burst 
position entered the LAT FoV again after To + 16.5 ks. 
No emission was detected and we derived upper limits 
for these late epochs as well. The Swift /XRT light curve 
is also plotted in the same figure. No flaring behavior is 
present in the X-ray light curve. 

The spectral properties of the GeV emission for differ- 
ent times are summarized in Table [TJ There is no signifi- 
cant spectral change between periods (I) a nd (HI) , as op- 
posed to preliminary results presented by iGruber et al.l 
(|2011f ). who found that the spectrum for the second emis- 
sion period was much softer than the first period (we were 
able to reproduce the same results by using Pass6 data). 
Since Pass7 data are recommended after their release, we 
will take the more conservative proposition that there is 
no subtle change in the spectral index in the LAT data. 

Figure [3] shows the spectral energy distribution of the 
X-ray and 7-ray emission integrated over the time inter- 
vals T + 255 s to T + 580 s and T + 220 s to T + 580 s, 
respectively. 

4. DISCUSSION 

The redshift of GRB 110625 A is unknown and we as- 
sume a redshift z ~ 1. As for the afterglow, the X- 
ray spectral index ~ —0.44 ± 0.22 and the early decline 
£-1.14 £ Q a tj me ~ 2 x 10 4 s) can be reproduced 
if the forward shock is in the slow cooling phase, i.e., 
Vm < < v c (where v m is the typical synchrotron 
radiation frequency and v c is the cooling frequency), 
the power-law index of the shock-accelerated electrons 
p ~ 2.5 and the circum-burst medium has a constant 
number density. In such a scenario, the X-ray data sug- 
gest that v c > 10 18 (i/2 x 10 4 s)" 1 / 2 . In the synchrotron 
radiation model, the cooling frequency is related to the 
cooling Lorentz factor 7 C as u c ~ m c cj^TB /2tt(1 + z), 
where m e is the rest mass of electrons, c is the speed 
of light, T » 88££,54ra" 1/8 (*/300 s )- 3 / 8 [(l + z)/2] 3 / 8 is 
the bulk Lorentz factor of the decelerating outflow and 
B k 3.5 Gauss e^E^n^it/^OO s)- 3 / 8 [(l + z)/2] 3 / 8 
is the strength of magnetic field in the emitting re- 
gion (e.g., Piran 1999). Please note that E^ is the ki- 



netic energy of the outflow, n is the number density of 
the circum-burst medium and £b is the fraction of the 
shock energy given to the magnetic field. The conven- 
tion Q n — Q/10 n has been adopted here and through- 
out this work. At t ~ 300 s, the time of the strong 
GeV emission, we have c c >8x 10 18 Hz and the corre- 
sponding cooling Lorentz factor 7 C ~ 2 x 10 6 (TB)~ 1,>2 ~ 

lO^B-^k 54 8n_1 ^ 8 - To account for the GeV spectrum 
one needs a typical energy of seed photons as low as 
- 100 MeV/7 2 - 2.4 x 10 12 Hz, which rules out the 
possibility that the GeV afterglow is dominated by the 
synchrotron self-Compton (SSC) radiation of the exter- 
nal forward shock. The synchrotron radiation of the 
forward shock electrons may produce prominent GeV 
emiss ion (e.g. JZou et al.l [20091 : iKumar fc Barniol Duranl 
12009ft . As shown in Fig.3, the GeV emission could be the 
high-energy tail of the synchrotron radiation component 
as long as the spectrum F v oc zy~(p _1 )/ 2 ^ ^-0.75 ^q^s 
in the energy range from 0.2 keV to ~ 100 MeV. Such 
a peculiar spectrum suggests a v c ~ 100 MeV, requiring 

Tvk~54 2e B 3 -2 n_1 ~ 3 x 10 6 , where we have taken the ex- 
pression of v c by Yost et al. (2003). Alternatively, one 
can imagine that there was a far-infrared flare and the 
forward shock electrons were cooled by the flare photons, 
producing a GeV radiation component via the inverse 
Compton scattering. These two kinds of processes may 
be able to account for some properties of the GeV emis- 
sion (e.g., the spectrum and the flux) but possibly not 
the temporal behavior. As shown in Fig.l, the possible 
abrupt decline of the GeV emission may be hard to be 
interpreted within the forward shock scenario. 

In the following investigation, we assume that the for- 
ward shock synchrotron GeV radiation is unimportant 

(requiring that v c <C 100 MeV, i.e., E^\j 2 e^^rT 1 <C 

3 x 10 6 ) and then examine the possibility that the GeV 
emission is the SSC radiation of an underlying ultra- 
violet (UV) flare with a peak energy -E Pj fl ar c ~ 20 
cV and a spectrum oc v~ x l for hv > T Pi fl arc (as al- 
ready mentioned in Sect. 3, there was no flaring be- 
havior in X-ray band). In the GRB afterglow, op- 
tical/ultraviolet flares could be powered via the so- 
called late int ernal shocks, as ob s erved for e xample in 
GRB 080129 (|Greiner et al.l[200l [G^o1[2009l) . Adopt- 
ing eq.(49) of Fan & Piran (2008), we have E s p sc ~ 

lOOMeV £:l /4 T fla 1 r / o 2 49 i? flarOil5 .5(T PiU v/20 eV) 2 , match- 
ing the data, where e = eB,flareAe, flare (where e e is the 
fraction of shock energy given to electrons in the flare 
phase), and the luminosity of the UV flare is related to 
the luminosity of GeV emission as Tfl arG ~ e 1 / 2 TccV ~ 
3 x lO 48 ^ 2 erg/s (where T GoV - 3 x 10 50 erg/s is 
the luminosity of the GeV emission). The requirement 
£ <C 1 suggests that the outflow is baryonic and the 
magnetic field in the emitting region is shock-generated. 
The required i?fl aro ~ 3 x 10 15 cm imposes a tight con- 
straint on the bulk Lorentz factor of the flare outflow, 
i.e., r flarc > [(l + z)i? flaro /160s/3x 10 10 cm s" 1 ] 1 / 2 ~ 25 
(The current data do not allow us to perform a reliable 
estimate on the variability timescale of the GeV emis- 
sion; that is why we used the whole period (see Fig.l) to 
constrain Ffl arc ). To be a valid interpretation, three more 
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Figure 1. Top panel: Weighted photon flux between 100 MeV and 20 GeV from GRB 110625A as observed using Fermi/LAT within a 
circular region of radius 10°, centered at the XRT position, between To + 180 s and To + 580 s. Each bin represents 40 s. The solid line 
represents the best-fit power law, with x 2 /d-o.f.= 16.79/8. The dashed line indicates the estimated number of background events averaged 
in the whole period. Bottom panel: Ratio, R := /~grb /^BG > °f the contribution from GRB 110625A over that from background based on 
a likelihood test for each time bin. The dotted line stands for R = 1. 



requirements should be satisfied: (i) the flare emission in 
X-ray band should be lower than the forward shock X- 
ray emission; (ii) 7 C > 10 5 still holds in the flare phase, in 
which the forward shock electrons suffer additional cool- 
ing by the UV photons; (iii) the GeV photons are not ab- 
sorbed via pair production on the high energy tail of the 
flare or on the SSC MeV photons. The first requirement 
is satisfied since the observed X-ray luminosity at t ~ 300 
s is ~ 3 x 1 48 erg/s while the emission of the UV flare in 
the X-ray band is - L flare (0.3 keV/£ p , flare )- - 7 ~ fO 48 
erg/s. The second requirement is met as long as the in- 
terstellar medium surrounding the burst has a number 
density n < O.Of cm~ 3 and ee < 0.01 (It is straight- 
forward to draw this conclusion with eq.(2) and eq.(5) 
of Fan & Piran 2006). Substituting the current physical 
parameters into eq.(43) of Fan et al. (2008), it is straight- 
forward to show that the absorption by pair production 
on the high energy tail of the flare can be ignored. The 
absorption by pair production on the SSC MeV photons 
can be calculated in a rather similar way and is found to 
be not serious 0. In other words, the third requirement 

2 With the SSC spectrum F E oc E' 1 / 2 for E < E* sc and F E oc 
E~ ( r T — for E > Ep SC , we just need to replace the expression of 
the term A r >E a in eq.(43) of Fan et al. (2008) by the new form 
87r(r 7 - 2)D£(£a£; sc )- 1/2 Fssc<^/[( 1 +^)(2r 7 - 3)] to estimate 
the pair production optical depth t 7 , where Fggc is the flux of the 



is also satisfied. Therefore we conclude that the SSC ra- 
diation of an underlying bright UV flare could account 
for the GeV radiation from t ~ 260 s to 340 s. The weak 
GeV emission in the time interval 460 — 500 s may be 
interpreted in the same scenario. 

5. CONCLUSIONS 

GRB 110625A is the third GRB detected by 
Fermi/LAT and Swift/XRT simultaneously. We have 
shown that the GeV light curve differs from a simple 
power-law decay, and probably consists of two emission 
periods. The rapid decrease of GeV flux during both 
periods challenges the notion that the emission comes 
from the external forward shock. While in the case of 
GRB 100728A, late-time X-ray flares seem to accompany 
the GeV emission, no such flares are seen in the time 
frame during which GeV emission was detected. This 
suggests a different origin of the GeV emission between 
the two cases. We discuss the possibility that the GeV 
emission is the SSC radiation of an underlying ultraviolet 

SSC radiation component and St is the observed typical variability 
timescale, _E a « 0.5 MeV rf 5 (B 7 /1 GeV) -1 is the energy of the 
photons that can absorb the GeV emission by pair production. 
For E^SC ~ 100 MeV, £ a - 0.5 MeV, i? flaro > 3 x 10 15 cm, 
F S sc ~ 10~ 7 erg s" 1 cm -2 , St ~ 80 s, T 7 ~ 2.7 and z ~ 1, we 
have t 77 Ri lla- T A r >E a /7207ri?| arc ~ 0.4 (e.g. Svensson 1987), 
where ctt is the Thompson cross section. 
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Figure 2. The LAT energy flu x derived from unbin ned likelihood analyses of GRB 110625A is shown in black. The XRT energy flux is 
shown in blue for comparison (Evan ct al. 2007, 2009). The X-ray light curve does not show any prominent X-ray flaring activity. 




Photon energy (eV) 



Figure 3. Spectral energy distribution of the X-ray and 7-ray emission integrated over the time intervals To + 255 s to To + 580 s and 
To + 220 s to To + 580 s, respectively. The butterfly represents the model absorbed power-law spectrum in the 0.3—10 keV range. 
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(UV) flare. Multiwavelength coverage of the rare class of 
LAT GRBs during the GeV-emitting period is crucial. 
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